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Abstract 

The self-consistent random phase approximation (RPA) approach with the residual interac- 
tion derived from a relativistic point-coupling energy functional is applied to evaluate the isospin 
symmetry-breaking corrections 5 C for the + — > + superallowed Fermi transitions. With these 5 C 
values, together with the available experimental ft values and the improved radiative corrections, 
the unitarity of the Cabibbo-Kobayashi-Maskawa (CKM) matrix is examined. Even with the con- 
sideration of uncertainty, the sum of squared top-row elements has been shown to deviate from the 
unitarity condition by 0.1% for all the employed relativistic energy functionals. 
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I. INTRODUCTION 

In the standard model, the Cabibbo-Kobayashi-Maskawa (CKM) matrix [l, 2] relates the 
quark eigenstates of the weak interaction with the quark mass eigenstates, and therefore it 
should be unitary. Examination of the unitarity of the CKM matrix provides a rigorous test 
for the standard model description of electroweak interactions. Usually, much attention has 
been paid to the top row of the CKM matrix. The unitarity requires, J2f=dsb l^«/| 2 = 1> 
where the values of \V US \ and \V u b\ can be taken from the Particle Data Group reviews js]. 
The value of matrix element \V u d\ can be determined in nuclear physics through the following 



ways: nuclear + — > + superallowed Fermi transition |4j, neutron decay [5|, pion f3 decay [6| 
and nuclear mirror transition jj. Among them, the first method provides the most precise 
determination of \V u d\ 31 5 given by \Gv/Gp\, where Gf is the Fermi coupling constant 
for purely leptonic decays. The vector coupling constant of semileptonic weak interactions 
Gy is determined by the nucleus-independent Tt value and transition-independent part of 
radiative corrections. Therefore, many efforts are devoted into the study of the nucleus- 
independent Tt value. 

On the theoretical side, there are several nuclear structure models that have been 
adopted to calculate the nucleus- independent Tt value, including the shell model Js | and 
the self-consistent charge-exchange random phase approximation (RPA) based on both non- 



relativistic [9] and relativistic |lOl . Illj energy density functionals, as well as the isospin- 
projection scheme on top of the Skyrme-DFT approach [3]. It has been shown in Ref. 11 | 
that the constancy of the Tt values can be satisfied for all the employed relativistic meson- 
exchange energy functionals. 

Recently, as the counterpart of meson-exchange energy functional, the relativistic point- 
coupling energy functional has attracted more and more attention because of its simplic- 
ity and success. In the point-coupling functional, there is no mesonic degree of freedom 
that makes its practical applications more feasible and the numerical effort considerably 
reduced [13|. In particular, the point-coupling functional can be easily extended to study 
the effects beyond the mean-field approximation that are important for a proper description 
of the nuclear low- lying collective excited states 14j-|l6l|. 

In view of these facts, recently, the model of change-exchange RPA with the residual 
interaction derived from the relativistic density-dependent point-coupling functional has 
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been developed [17|. In this work, we will adopt the similar model, but we'll start from 
the non-linear version of point-coupling functional to study the isospin symmetry-breaking 
corrections 8 C in the + — » + superallowed transitions. With these 5 C values, together with 
the most recent experimental ft values 18] and the improved radiative corrections js), I19I ]. 
the unitarity of the CKM matrix will be examined. 

This paper is arranged as follows. The model is introduced briefly in Sec. [Til The 
results and corresponding discussions are given in Sec. IIII1 A brief summary of the present 
investigation is presented in Sec. |TVJ 

II. THE MODEL 

For the charge-exchange channels of both the r + and t__ , the relativistic RPA equation 
has the following form [ll]. 



A 3 - E J - 

t ^pn,p'n pn,n'p' 

—B J - — A J - 

np,p'n' '^ l np,n'p l 















\ np 



(1) 



In the above equation, p and p {n and n) denote occupied and unoccupied proton (neutron) 
states, where the unoccupied states include the states above the Fermi surface and those in 
the Dirac sea. U Ju and V Jv are the amplitudes corresponding to the RPA energy oj u . The 
matrices A and B are given by, 

^12,34 = (£1 - £2)^12,34 + (14|V ph |32), (2) 
£12,34 = -(13|y ph |42), (3) 

where Ei,E 2 are the single-particle energies from mean-field calculations. 

For the sake of the self-consistency, the particle-hole residual interaction Vph is de- 
rived from the same effective Lagrangian density as the mean-field Dirac single-nucleon 
Hamiltonian that determines the mean-field for the nuclear ground-state. The spin-isospin- 
dependent interaction terms are generated by the isovector terms. Although the direct 
one-7r contribution to the nuclear ground-state vanishes at the Hartree level because of par- 
ity conservation, the pion nevertheless must be included in the calculation of spin-isospin 
excitations 20]. 
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The derivative type of the pion-coupling necessitates the inclusion of the zero-range 
Landau-Migdal term, which accounts for the contact part of the nucleon-nucleon interaction, 

V s ^(l, 2) = g'[— 7o757^*]i[— lolhl^iri - r 2 ), (4) 

where /^/47r = 0.08, m n = 138 MeV. The parameter g' in principle has to be adjusted to 
reproduce the experimental Gamow- Teller resonance excitation energy |20|. However, the 
direct contribution from the 7r-meson field vanishes for + — > + transition. Therefore, the 
self-consistency is still maintained in the RH+RPA study of superallowed Fermi transition 



and we take g' = 1/3 as Ref. U| 



In the captions, the re.ativistic point-coupling energy fnnetiona. of both PC Fl M 
and PC-PK1 [13J is adopted, where the scalar-isovector coupling terms do not exist. As a re- 
sult, the particle- hole residual interaction V^h that gives non-zero contribution to the charge- 



exchange RPA matrix is composed of three parts: the vector-isovector coupling terms 



22|, 



V TV (1, 2) = [ 7 o7^i(«w + <W A) [707^25(17 - r 2 ), (5) 
the Landau-Migdal term (TJJ, and n- meson field term fill. I^. 

K(l,2) = -[— 7q757 ■ Vf]i[— 70757 ■ Vf] 2 Ar(ri, r 2 ), (6) 
where Z) 7r (r 1 ,r 2 ) denotes the 7r-meson propagator, 

1 p-m w |ri-r 2 | 

otTv an d Stv are the coupling; constants in the vector-isovector channel of point-coupling 
energy functional 
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21] , and A is the Laplace operator. 7 M are the four- component Dirac 
matrices. At the Hartree level, the Coulomb term does not contribute to the charge-exchange 
RPA matrix element. 

The eigenvectors of the RPA Eq. (CQ) are separated into two groups, which respectively 
represent the excitations of the r_ and r + channels with the following normalization condi- 
tions 

W ) 2 - E(^t) 2 = +1 > for r ~ channel > ( g ) 

pn np 

E^) 2 - E(^) 2 = - 1 ' for r + channel • ( 9 ) 



pn np 
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The excitation energies Q u and the corresponding forward X Ju , backward Y Jv amplitudes 
in the r_ and r + channels are given by, 



Vl u = —ui u , X, 



vJv 

pn 


= U J - 

^ pn ' 


yju 

np 


= V J - , 

' np ' 


vJv 
np 


= V J -, 

np ' 


Y Jv 

pn 


= U J - 

pn 



for r_ channel 



for r + channel 



(10) 
(11) 



Subsequently, it is straightforward to calculate the + — >■ + superallowed transition proba- 
bilities between the ground-state and excited states with the forward and backward ampli- 
tudes for the r_ and r+ channels as follows, 

2 

,, _\ . ^ ^ /_ n n v 

r_ n 



B 



Jv 



B 



with k 



ab 



Jv 



\3a+3b 



J2 X in(P II T- || n) +J2 K nP Y nP<S> II r- || n) 

pn, np 

J2 X np(n II r+ || p) + ^KpnY^in || r + \\ p) 



np 



pn 



(12) 
(13) 



III. RESULTS AND DISCUSSIONS 

For simplicity, the pairing correlations are neglected and the filling approximation is 
used. Moreover, the calculation is presently restricted to preserve spherical symmetry. As 
a result, the Dirac equation for nucleons can be solved easily in coordinate space using the 



numerical techniques in Ref. 23|], where the box size R = 15 fm and the mesh size Ar = 0.1 
fm. The solutions of Dirac equation, including single-particle wave functions and energies, 
are used as inputs of the RPA equation. As usual, the single-particle energy truncation is 
introduced. We find that with the choice of the truncations [—940, 1100] MeV for the PC-F1 
and [—940, 1160] MeV for the PC-PK1, the model- independent sum rule 



N — Z 



(14) 



can be fulfilled up to 10~ 5 accuracy, and the isospin symmetry-breaking corrections 5 C are 
stable with respect to these numerical inputs at the same level of accuracy. 

The isospin symmetry-breaking correction S c is determined by the superallowed transition 
strength Mp as follows, 

M|^|(/|r ± K)| 2 = |M | 2 (l-5 c ), (15) 
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FIG. 1: Distribution of the strengths B u (0 + — > 0+) for the superallowed transition °C — >• B. 

where Mo = \/2 for T = 1 states with the exact isospin symmetry, and |/), \i) represent 
the final and initial states of t± transitions respectively. It has to be pointed out that in 
practical calculations, there are several final states 0+ with different excitation energies E v 
and strengths B v as shown in Figure [1] for the superallowed transition, for example, 10 C 
— > 10 B. As usual, the final state |/) is taken to be the excited state with the largest transi- 
tion strength B v (close to 2). Table [J presents the isospin symmetry-breaking correction S c 
in several 0+ — > 0j superallowed transitions, calculated using the RPA for the relativistic 
point- c oup ling energy functional of both PC-F1 and PC-PK1. For comparison, the results 
in Ref . [ll] from the RPA calculations but with the relativistic meson-exchange energy func- 



tionals PKOl 24] and DD-ME2 25] are given as well. It has been found in Ref. [llj that the 
significant differences in 5 C values obtained by PKOl and by DD-ME2 are mainly due to the 
different treatments of the Coulomb field. The Fock terms have been included in the former 
case (PKOl), but not in the latter case (DD-ME2). Table [J shows that the point-coupling 
functionals PC-F1 and PC-PK1 without the Fock terms give quite similar S c values, both of 
which are in between the results by the PKOl and DD-ME2. 

In Table [Til the excitation energies E v for the Of — > 0j superallowed transitions are 
shown. For comparison with the experimental values taken from the recent survey 18] , 
the corrections due to the proton-neutron mass difference in particle-hole configurations are 
made for the calculated results. It is shown that all the relativistic energy functionals can 
reproduce the excitation energies quite well. 
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TABLE I: Isospin symmetry-breaking corrections S c (in %) for the + — > + superallowed transi- 
tions obtained by RPA with the relativistic point-coupling energy functional of both PC-F1 and 
PC-PK1. The results are compared with those obtained with the relativistic meson-exchange 
energy functionals PKOl and DD-ME2 |n|. 



PC-F1 PC-PK1 PKOl [11] DD-ME2 [llj 



io c _ 


^ 10 B 


0.109 




0.082 


0.150 


14Q- 


^ 14 N 


0.150 




0.114 


0.197 


18 Ne- 


^ 18 F 


0.309 


0.297 


0.270 


0.430 


26 Si- 


^ 26 A1 


0.202 


0.180 


0.176 


0.252 


30 S- 


^30p 


0.420 


0.488 


0.497 


0.633 


34 Ar _ 


^ 34 C1 


0.379 


0.378 


0.268 


0.376 


38 Ca- 


^38 K 


0.347 


0.325 


0.313 


0.441 


42 Ti _ 


* 42 Sc 


0.400 


0.375 


0.384 


0.523 


26 Al- 


+ 26 Mg 


0.159 


0.141 


0.139 


0.198 


34 C1 _ 


^ 34 S 


0.316 


0.314 


0.234 


0.307 


38 K _ 


-)■ 38 Ar 


0.294 


0.275 


0.278 


0.371 


42 Sc- 


^ 42 Ca 


0.345 


0.322 


0.333 


0.448 


54 Co- 


^ 54 Fe 


0.339 


0.301 


0.319 


0.393 


66 As- 


^ 66 Ge 


0.522 


0.488 


0.475 


0.572 


70 Br _ 


^ 70 Se 


0.935 


0.998 


1.140 


1.268 


74 Rb- 


^ 74 Kr 


0.668 


0.597 


1.088 


1.258 



The nucleus-independent Tt can be obtained by the experimental ft value, isospin 

symmetry-breaking correction 5 C and transition-independent part of radioactive correc- 
tions |4J, 

Ft = ft(l + 5' R )(l + 5 NS -5 c ), (16) 

where 5' R is the part of nucleus-independent radiative correction, which is a functional only of 
the electron's energy and the charge of daughter nucleus Z while 5ns is the part of radiative 
correction that depends on the details of nuclear structure. 



With the recent theoretical results of 5 R and 5ns 



and the recent experimental ft 



TABLE II: Excitation energies (in MeV) for the + —> + superallowed transitions measured by 
taking the ground-state of the corresponding even-even nuclei as reference. For comparison with 
the experimental values taken from the recent survey [3] , the corrections due to the proton- neutron 
mass difference in particle-hole configurations are made for the calculated results. 







Expt. 


PC-F1 


PC-PK1 


PKOl [111 


DD-ME2 fill 


10 n - 

V 


4 10 R 




9 917 






- . LixJXJ 


14 - 


■> 14 N 


-2.831 


-2.967 




-2.420 


-3.081 


18 Ne- 


■> 18 F 


-3.402 


-3.400 


-3.419 


-3.195 


-3.451 


26 Si- 


* 26 Al 


-4.842 


-4.956 


-5.055 


-4.531 


-5.110 


30 S - 


v 30p 


-5.460 


-5.295 


-5.330 


-4.845 


-5.395 


34 Ar - 


■»• 34 C1 


-6.063 


-5.975 


-5.964 


-5.559 


-6.278 


38 Ca- 


■»• 38 K 


-6.612 


-6.582 


-6.576 


-6.035 


-6.775 


42 Ti- 


+ 42 Sc 


-7.000 


-6.833 


-6.869 


-6.661 


-6.964 


26 A1 _ 


■»• 26 Mg 


4.233 


4.211 


4.304 


3.908 


4.350 


34 C1 _ 


^ 34 S 


5.492 


5.292 


5.288 


5.062 


5.561 


38 K _ 


* 38 Ar 


6.044 


5.905 


5.905 


5.557 


6.083 


42 Sc- 


^ 42 Ca 


6.426 


6.207 


6.242 


6.118 


6.333 


54 Co- 


^ 54 Fe 


8.244 


8.016 


8.122 


7.720 


8.240 


66 As- 


^ 66 Ge 


9.579 


9.367 


9.434 


9.044 


9.677 


70 Br _ 


^ 70 Se 


9.970 


9.735 


9.806 


9.632 


9.852 


74 Rb- 


^ 74 Kr 


10.417 


10.246 


10.299 


10.005 


10.437 



values [18], the nucleus- independent Tti value and its uncertainty bTt% for each superallowed 
Fermi transition are obtained with the calculated 8 C values and listed in Table 1111} where 
the corresponding average value J-t and \ 2 l v are given as well. In the calculation of &Tti, 
the uncertainty of 8 C is taken as zero. The average value Tt is calculated in the following 
way |26(, 

Ti±5Ti= I^Ill ± (V c^)- 1/2 (17) 
l^i^i - 

where the weight is 

ut = l/^FUf. (18) 



8 



The x 2 l v is calculated by 

- Tf - Ft 

^ = E(^f ! ) 2 /(^-i) (19) 

i 

where N is the number of the calculated superallowed Fermi transitions. Since the un- 
certainties 5Fti in the superallowed Fermi transitions 34 Ar — > 34 C1 and 74 Rb — > 74 Kr are 
obviously larger than those in other transitions and excluded in the calculations of average 
value Ft and x 2 l v - 

TABLE III: Nucleus-independent Ft values (in s), its average value Ft (in s) and the x 2 l v - 







PC-F1 


PC-PK1 




io B 


3078.7(45) 




14 0^ 


14 N 


3077.0(31) 




34 Ar ^ 


34 C1 


3078.5(84) 


3078.5(84) 


26 Al^ 


26 Mg 


3077.0(13) 


3077.6(13) 


34 C1^ 


34 S 


3081.0(15) 


3081.1(15) 


38 K ^ 


38 Ar 


3083.6(16) 


3084.2(16) 


42 Sc^ 


42 Ca 


3082.3(21) 


3083.0(21) 


54 Co^ 


54 Fe 


3083.3(24) 


3084.4(24) 


74 Rb^ 


74 Kr 


3119.6(88) 


3120.2(88) 


averag 




3080.3(7) 


3081.1(7) 


X 2 /v 


1.1 


1.4 



Table HIl shows that the \ 2 jv is 1.1 for PC-F1 and 1.4 for PC-PK1. It indicates that the 
constancy of the Ft values is good for both cases. 

To illustrate the constancy of the Ft values more clearly, we plot the nucleus-independent 
Ft values as a function of the charge of the daughter nucleus Z in Figure [2J The shaded 
horizontal band gives one standard deviation around the average Ft value. 

With the nucleus-independent Ft value, the absolute value of element V uc i in the CKM 
matrix can be calculated by 

K 

|Kd|2 = 2Gl{l + Al)Fi m 
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FIG. 2: Nucleus-independent J-t values as a function of the charge of the daughter nucleus Z. The 
shaded horizontal band gives one standard deviation around the average Ft value. 

where K/(Kcf = 8120.2787(11) x 10~ 10 GeV -4 s, the purely leptonic decays Fermi coupling 
constant GV/(/ic) 3 = 1.16637(1) x 10~ 5 GeV -2 [3[] and the radiative corrections' transition- 



independent part A£ = 2.361(38)% 

Together with the other two CKM matrix elements \V US \ = 0.2255(19) and \V u b\ = 
0.00393(36) [3], the unitarity of the CKM matrix can be examined by the sum of squared 
top- row elements of the CKM matrix. The values of \V U d\ and |V U( i| 2 + |Kts| 2 + |Ktb| 2 are listed 
in Table IIV] where we also give the values of uncertainty, which might be underestimated 
to some extent as the uncertainty of 5 C was assumed to be zero and the systematic errors 
were not taken into account. 

TABLE IV: The absolute value of matrix element V u d and the sum of squared top-row elements of 
the CKM matrix from the RPA calculations with different relativistic energy functionals. 



\V ud \ 



\v ud \ 2 + \v us \ 2 + \v ub {< 



PC-F1 
PC-PK1 
PKOl M 



0.97290(21) 
0.97278(22) 
0.97273(27) 



DD-ME2 [11] 0.97311(26) 



0.9974(10) 
0.9972(10) 
0.9971(10) 
0.9978(10) 



Table HV] shows that the \V u d\ values obtained by PC-F1 and PC-PK1 are quite similar 
and close to the results of PKOl. In addition, it is seen that even with the uncertainty, 
the sum of squared top-row elements deviates from the unitarity condition by 0.1% in the 
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charge-exchange RPA calculations with these four relativistic effective interactions. 



IV. SUMMARY 



In summary, the self-consistent RPA with the residual interaction derived from the rela- 
tivistic point-coupling energy functional has been applied to calculate the isospin symmetry- 
breaking corrections 5 C for several typical + — > + superallowed transitions. Together with 
the experimental ft values in the most recent survey and the improved radiative corrections, 
the corresponding nucleus- independent J-t values and matrix element jV^I have been cal- 
culated. It has been found that the \V u d\ values obtained by PC-F1 and PC-PK1 are quite 
similar and close to the results of PKOl. However, even with the uncertainty, the sum of 
squared top-row elements has been shown to deviate from the unitarity condition by 0.1% 
for all the employed relativistic energy functionals. It indicates that other effects, includ- 
ing deformation and pairing correlations, would play important roles. Therefore, it is very 
interesting to study the + — » + superallowed transitions in the framework of deformed 
QRPA with a proper pairing force, for instance a separable pairing force 27|. Of course, in 
this case, particle number projection and angular momentum projection are required to give 
good nucleon number and angular momentum for RPA states. Work along this direction is 
in progress. 
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